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Abstract 

The buckling loads of laminated plates are predicted using a new theory which 
takes into account transverse shearing effects. This new theory assumes 
trigonometric terms through-the-thickness in the displacements to take into 
account transverse shearing effects in thick plates. Buckling loads predicted by 
the new theory and by traditional theories are compared for isotropic and 
laminated plates. The effect of ply orientation on the buckling loads predicted 
by each theory is demonstrated. 

Introduction 

Classical plate theory accurately predicts the inplane stresses, displacements 
and buckling loads of thin plates of homogeneous materials. Conventional 
transverse shear deformation theory accurately predicts the stresses, 
displacements and buckling loads of thicker plates of homogeneous materials. 
However, these theories are not as accurate when applied to laminated plates as 
they. are when applied to isotropic plates because they do not adequately take 
into account transverse shearing effects. Since transverse shearing can 
significantly affect the response of laminated plates, a theory which predicts 
the buckling loads and transverse shearing effects more accurately than 
traditional theories is essential for predicting the behavior of laminated 
plates. Such a theory is developed in reference 1. 
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The purpose of the present paper is to apply the theory developed in reference 1 
to simply supported flat plates and to compare the buckling loads for the theory 
of reference 1 , for classical theory and for the conventional transverse shear 
deformation theory for plates of different thicknesses. In addition, the effect 
of changes in ply orientation on buckling load for thick laminated plates is 
discussed. 

Analysis Approach 

The results presented in this paper are obtained by applying the theory developed 
in reference 1 * In reference 1 the potential energy method is used to obtain 
equations for the buckling of a flat plate with the coordinate system and plate 
dimensions shown in the sketch. 



The series chosen to represent the plate displacements are shown in equations (1) 
and have three types of terms. 

e 

u(x,y,z) * [u° + u a (z/h) + u^ sin(irz/h) ] cos (mix /a) sin(miry/b) 
v(x,y,z) - [v° + v a (z/h) + v 1 sin(irz/h)] sin(mrx/a) cos(miry/b) (1) 

w(x,y,z) * [w° +W 1 cos(itz/h)] sln(nirx/a) sin(miry/b) 
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The traditional terms from classical plate theory (those independent of position 
in the z direction), and from conventional transverse shear deformation theory 
(those linear in z) are represented by the terms with superscripts o and a, 
respectively. The assumed displacement series also includes trigonometric terms 
in the through-the-thickness or z direction, which are represented by the terms 
with superscript 1. By including trigonometric terms in the assumed displacement 
series, a more accurate solution is obtained because more transverse shear 
deformation is permitted than in the solutions with fewer assumed displacement 
terms. 


Differential equations for predicting the buckling loads of flat plates are 
developed in reference 1 by using the potential energy method and the assumed 
displacement series shown above. The equilibrium equations for the flat plate 
are obtained from: 
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Per* the ease Of w ■ 0 before buckling* the differential equations of equilibrium 
can be converted into buckling equations by considering small changes in 
deformations that occur at buckling. This assumption reduces the third of 
equations (2) to: 


4/2 Ct xx,x * V.y 3 dz * loa<lin8 tePms ' 0 (3) 

while the other equations of equations (2) remain unchanged. 

The differential equations represented by equations (2) and (3) can be reduced to 
linear equations containing only the displacement and loading terms as unknowns 
by using the following stress-strain and strain-displacement relations 
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where [C^] is the material stiffness matrix. 


Inserting equations (4) into equations (2) and (3) and integrating in the 
thickness direction gives the following linear equations. 
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The displacements (u°, u a , u\...) in equations (5) are now the small changes 
which occur at buckling and n and m are the number of half waves in the x and y 
directions, respectively. 


2 2 o 

Loading terms in the form {N (mr/a) + N (mir/b) }w for compressive buckling are 

x y 

added to the third of equations (5). Buckling loads can be found from these 
equations by reducing equations (5) (with the loading terms) to the form of 
equation (6), 
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The lowest value of N or N for which the determinant of [K] - 0 is the buckling 

x y 

load. The following section contains numerical results for the compressive 
buckling of simply supported plates loaded in one direction. 

Results and Discussion 

A comparison is made between the results from classical plate theory in which no 
transverse shearing is assumed, from conventional transverse shear deformation 
theory in which the transverse shear stress is assumed to vary only along the 
length and width of the plate, and from the theory of reference 1 in which the 
transverse shear stress varies in all three coordinate directions. The 
conventional transverse shear deformation theory can be obtained by neglecting 
the trigonometric terms In the theory of reference 1 and the classical theory can 
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be obtained by assuming the plate to be very thin. These three theories are 
applied to an isotropic plate, a [*45]^ laminated plate and a [0].^ laminated 

plate. All of the results presented are for square plates. A variety of 
thicknesses ranging from very thin to one-third the length are examined. The 
buckling loads and stresses at buckling predicted by each theory are compared. 

The relationship between ply orientation and buckling load is also studied for 
[ ± 0] 4s laminates for several thicknesses and for ply orientations ranging from 0 

to 90 degrees. 

The buckling load parameter, k = a N / ■( *' /D^ , is shown in figure 1 for a 

square plate with thicknesses varying from 0 to .3 times the length for the three 
theories. In the case of the thin isotropic plate, the difference between the 
buckling loads predicted by each theory is less than one percent. The classical 
theory, which assumes no transverse shearing (t xz = T yz = 0), gives a buckling 

load parameter which is independent of thickness and therefore only accurate for 

thin plates. The theory in which the transverse shear stress is assumed to be 

constant through-the-thickness (t (x,y); t (x,y)) gives a buckling load which 

xz yz 

is slightly lower than the buckling load from classical theory for thin plates 

and a buckling load which is 30 percent lower for very thick plates (t/a = .3). 

The theory of reference 1 , in which the transverse shear stress is a function of 

all three coordinates (t (x,y,z) ; t (x,y,z)) , gives a buckling load which is 36 

xz yz 

percent lower than the buckling load from classical theory for the thickest plate 


studied 
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The buckling load parameter versus the thickness-to-length ratio is shown fop a 


laminated plate of [±45]^ ply orientation in figure 2, 


Here the same trend 


occurs as fop the isotropic piste. Fop very thin pistes all theories give 
essentially the same buckling loads. But for a plate which has a thickness-to- 
length ratio of .3. the convent ional tpans verse shear deformation theory and the 
theory of reference 1 have predicted buckling loads that are 22 and 36 percent 
lower, respectively, than results from classical theory. 


The buckling stresses versus the thickness-to-length ratio for a [+45]^ 

laminated plate and for a [0]^ plate are shown in figure 3. While the addition 

of the trigonometric terms decreases the predicted buckling stress in both cases, 
the more significant decrease is caused by including the trigonometric terms in 
the [±453^ laminate. 


A comparison of the buckling load parameter for the three theories is shown in 

figure 4 for ply orientations ranging from 0 to 90 degrees. The plate3 examined 

here have thickness-to-length ratios (t/a) of .05 for a very thin plate, .20 for 

a plate with intermediate thickness, and ,30 for a very thick plate. The 

classical theory gives the same curve as the conventional transverse shearing 

theory for the thinnest plate. Since no thickness effects are taken into account 

in the classical theory, all thicknesses give the same curve. There is almost no 

difference in the buckling loads predicted by each theory for the very thin plate 

for all ply orientations. The difference between the conventional transverse 

shearing theory (t (x,y) and t (x,y) are constant through-the-thickness) and 

xz yz 

the theory of reference 1 for the plate with intermediate thickness (t/a=.20) 
varies from 5 percent for plates with 0 = 45° to 13 percent for 0 = 30°. For the 
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thick plate (t/a = .30) the reduction in buckling load is greater than 12 percent 
for all orientations for the theory of reference 1 compared to the conventional 
transverse shear theory. The largest difference in buckling load is a reduction 
of 17 percent for a plate with 0 = 45°. In all cases except the very thin plate, 
the buckling loads predicted by the theory of reference 1 are lower than those 
predicted by the classical theory and the conventional transverse shear 
deformation theory. The theory of reference 1 allows more flexibility than the 
other two theories by including more terms in the assumed displacement series. 

It is a higher order theory and is, therefore, more accurate. 

Concluding Remarks 

Results are obtained for a new theory for predicting the buckling load of simply 
supported laminated plates. The new theory uses trigonometric terms in the 
assumed displacement series to predict the effects of transverse shearing. 
Transverse shearing causes significant decreases in the buckling load results of 
thick laminated plates. For thick plates the new theory predicts buckling loads 
which are as much as 30 percent lower than those predicted by traditional 
theories. 

A parametric study indicates that the ±45° ply orientation is the orientation 
most affected by transverse shearing in thick laminates. A comparison of the 
buckling load results from a theory which does not include the trigonometric 
terms (conventional transverse shear deformation theory) to the results from the 
new theory shows a decrease in buckling load results from 12 to 17 percent, 
depending on the ply orientation. 



Stein, M.s Nonlinear Theory for Laminated and Thick Plates and Shells 
includihg the Effects of Transverse Shearing, Presented at the 26th 
AIAA/ASME/ASCE/AHS Structures , Structural Dymani cs and Materials Conference, 
Orlando, Florida, April 15-17, 1985. AIAA paper 85-0671 -CP, April 1985. 
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Figure 1. Comparison of buckling load parameter of an 
isotropic square plate as a function of thickness—to— length 
ratio obtained by several aoDroximations . 
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Figure 2. Comparison of buckling load parameter of a 
laminated graphite-epoxy, [±45] 4g , square plate as a 

function of thickness-to-length ratio obtained by several 
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Figure 4. Comparison of buckling load parameter f 
graphite-epoxy laminated plates as a function of ply orie 
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